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ABSTRACT 

Aims. We study a new mechanism for deflagration to detonation transition in thermonuclear supernovae (SNe la), based 
on the formation of shocks by amplification of sound waves in the steep density gradients of white dwarfs envelopes. 
We characterize, in terms of wavelength and amplitude, the perturbations which will ignite a detonation after their 
amplification. 

Methods. This study was performed using the well tested HERACLES code, a conservative hydrodynamical code, 
validated in the present specific application by an analytical description of the propagation of sound waves in white 
dwarfs. Thermonuclear combustion of the carbon oxygen fuel was treated with the a-chain nuclear reactions network. 
Results. In planar geometry we found the critical parameter to be the height of shock formation. When it occurs in the 
inner dense regions (p > 10® g.cmT'^) detonation is inevitable but can take an arbitrarily long time. We found that 
ignition can be achieved for perturbation as low as Mach number: M~ 0. 005, with heating times compatible with typical 
explosion time scale (a few seconds). On the opposite no ignition occurs when shocks initiated by small amplitude or 
large wavelength form further away in less dense regions. We show finally that ignition is also achieved in a spherical 
self-gravitating spherical model of cold C-l-O white dwarf of 1.430 Mq, but due to the spherical damping of sound 
waves it necessitates stronger perturbation {M~ 0.02). Small perturbations {M~ 0.003) could still trigger detonation 
if a small helium layer is considered. In the context of SNe la, one has to consider further the initial expansion of the 
white dwarf, triggered by the defiagration, prior to the transition to detonation. As the star expands, gradients get 
flatter and ignition requires increasingly strong perturbations. 
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1. Introduction 

Type la supernovae are thought to result from the ther- 
monuclear explosion of a carbon oxygen white dwarf. 
Presently two main families of models are proposed: in the 
single degenerate scenario (Whelan & Iben 1973) explo- 
sion is triggered by continuous accretion of the hydrogen or 
helium rich envelop of a non-degenerate companion until 
Chandrasekhar mass is reached; in the double degenerate 
scenario (Iben & Tutukov 1984) the explosion results of 
the merging of two C+0 white dwarfs. In the latter mech- 
anism, only the detonation mode of combustion can release 
enough energy to prevent a collapse and produce a healthy 
explosion. Furthermore ignition is not so easy and require 
violent mergers (Pakmor et al. 2012), that might not be 
frequent enough to explain the SN la rates. In the former 
scenario pure detonation and pure deflagration fail to repro- 
duce light curves and nucleosynthesis. Detonation models 
(Arnett 1969) incinerate the whole star at high density pro- 
ducing mostly ^^Ni and resulting in wrong nucleosynthe- 
sis. On the other hand pure deflagration is very subsonic 
and gives time for expansion so that burning also occurs 
at lower density and produce the intermediate elements 
needed for spectra, but those models are under-energetics 
(Ropkc ct al. 2007). Accordingly the delayed detonation 
scenario has been proposed (Khokhlov 1991a) to combine 
the advantages of both combustion modes. In this model 



a flame born as a subsonic deflagration is accelerated by 
turbulence and transitions at some point to a detonation. 
The initial phase allow the star to expand so that IME are 
produced in the final detonation stage. 

However, after many years of studies the mechanism 
of deflagration-to-detonation transition (DDT) is still un- 
clear. It has been shown by 3D simulations (Gamczo ct al. 
2005) that this mechanism was the most promising path for 
single degenerate models to reasonably account for observ- 
ables, however due to the lack of theoretical understanding 
of this transition, the turning point between deflagration 
and detonation still relies on educated guesses (Woosley 
ct al. 2009) or is adjusted in order to best reproduces nu- 
cleosynthesis of normal type la supernovae (Jackson ct al. 
(2010)). On earth, DDT has been observed in confined me- 
dia, for which Ivanov ct al. (2010) provided numerical ev- 
idence that the transition may occur through shock-fiame 
interactions, see also Khokhlov et al. (1999) and Oran & 
Gamczo (2007) for more details on shock-flame interaction 
in DDT. See also CiccarcUi & Dorofccv (2008) for a review 
of DDT mechanisms in typical terrestrial systems. However 
in unconfined media, terrestrial experiments showed that a 
turbulent flame is rather stable against DDT. Thus it has 
been proposed that the transition should occur through 
the gradient mechanism (Khokhlov 1991b), where a suit- 
able gradient of the reactions induction time is created and 
maintained over a time long enough for a coherent burn- 
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ing wave to build up pressure and create a self-sustained 
detonation. It has been suggested that the turning point 
from turbulent deflagration to detonation might be the cre- 
ation of such a flat gradient at the transition between the 
flamelet and the distributed combustion regimes, for which 
the flame changes in nature. However it is not clear if the 
turbulence level necessary for the distributed regime to oc- 
cur is achieved during the explosion, even at the late stages. 
More recently, through numerical studies Poludnenko et al. 
(2011) claimed that intense enough turbulence could di- 
rectly initiate DDT, without relying on this induction time 
gradient. Other mechanisms able to ignite a delayed detona- 
tion from a turbulent deflagration have also been proposed: 
creation of hot spot, or shock convergence in the turbulent 
flame brush (Kushnir et al. 2012). Thus in the context of 
supernova where the plasma is fully unconfined but where 
high level of turbulence should be present, the question is 
still debated. 

Here we examine a new mechanism by which pressure 
perturbations, created by flame turbulence, are amplified 
during their propagation ahead of the flame in the steep 
density gradient where they turn into shocks, which can 
be strong enough to ignite a detonation. It is known since 
Lighthill (1955) that compressible turbulence can produce 
sound waves, even with small rms Mach number turbu- 
lence. Moreover the turbulent combustion is likely to en- 
hance sound production. It is also known that sound waves 
propagating through steep density gradients get amplified, 
and given a large enough density jump they result in strong 
shocks and substantial heating. The mechanism of forma- 
tion of strong shocks in density gradient has been proposed, 
for example, for the heating of the solar chromosphere (see 
Ulmschncidcr (1970) for a review). 

In this paper, we analyse the propagation properties of 
finite amplitude pressure waves in C-l-0 degenerate matter, 
first in terms of the density gradient and then in a more re- 
alistic spherical WD model. We determine the frequencies 
and amplitudes which lead to a successful detonation. We 
postpone to a following paper the study of the actual source 
of pressure perturbations. Our paper is constructed as fol- 
lows: in the first part we present our numerical methods 
and the physics included. In the second part we derive an 
approximate equation for shock formation and amplifica- 
tion in density gradients in C+0 degenerate matter and 
compare numerical results obtained for the planar models 
to validate the numerical simulations. Then we apply this 
to spherical models closer to the reality: first a cold and 
dense white dwarf and then a set of pre-expanded struc- 
tures. Finally we briefly discuss the applicability of this 
mechanism to thermonuclear supernovae. 



2. Numerical Methods and Physical Input 

We will consider the dense plasma of WD interiors com- 
posed of carbon - oxygen mixture. At these densities 
the plasma is degenerate and needs special treatment. 
Subsection 2.1 describes the method used to simulate gas 
hydrodynamic in the dense plasma with the HERACLES 
code^ (Gonzalez et al. 2007). Subsection 2.2 describes the 
implementation of the degenerate equation of state (EoS) 



within the code. Finally subsection 2.4 presents the imple- 
mentation of the thermonuclear reaction network. 



2.1. The Code: HERACLES 

The simulations are performed using the HERACLES code. 
It is a grid-based conservative compressible hydro code, 
with a second order Godunov scheme, which solves the 
usual Euler fluid equations: 
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where g is the acceleration of gravity, enuc is the energy 
release due to thermonuclear reactions and E = Eo + Etii + 
0.5pu^ is the total energy including respectively degener- 
acy, thermal and kinetic energy. The fluxes between cells are 
computed through a Riemann solver adapted to our degen- 
erate EoS. It also include gravity implemented as a Cranck- 
Nicholson step decoupled from the hydro step. Since we 
made only ID calculations, the gravitational acceleration is 
either constant in the planar geometry or is computed self- 
consistently at every radius r as the acceleration of gravity 
due to the enclosed mass: 

Gm(r) 
9i^> = 

where m{r) — Airr^ p{r)dr is the enclosed mass in the 
sphere of radius r and G is gravitational constant. 

2.2. Equation of State 

In a white dwarf, the electron gas is degenerate and ions 
form an almost perfect gas. Depending on the tempera- 
ture and density, electrons have various degrees of degen- 
eracy. Thus a general EoS is used to described the plasma 
including a perfect gas of ions, photons in local thermal 
equilibrium and electrons at all degrees of degeneracy and 
relativity (Timmes & Swcsty 2000). To solve the Riemann 
problem for a general EoS we use the solver described by 
Colella & Glaz (1985) as implemented in the FLASH code 
(Fryxell et al. 2000). But in some cases it was not accu- 
rate enough to follow precisely the temperature. Indeed 
in the density and temperature conditions relevant here, 
the main contribution to the pressure and energy of the 
gas is from the degenerate electrons and is essentially in- 
dependent of the temperature. The thermal contribution, 
defined as Eth{p,T) = E,r,{p,T) - E,„{p,T = OK) rep- 
resents often less than 0.1% of the internal energy [Em). 
Since a conservative code computes only the total energy: 
E = Eo + Eth + 0.5pit^, direct inversion to obtain the tem- 
perature is sometimes inaccurate. Thus, in order to accu- 
rately follow the temperature of the plasma, an equation of 
evolution for the thermal energy is derived and used when 
necessary: 
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where Pth is the thermal pressure defined as Pth {p, T) = 
Pip,T)-P{p,T^OK). 

2.3. Numerical Dissipation and Waves Propagation 

In order to resolve sound waves with minimum numerical 
dissipation, there must be enough points per wavelength 
and the time step must be much smaller than the period. 
To check the ability of our code to propagate sound waves 
without too much dissipation, we initiate waves of different 
frequencies at a given amplitude and measure the ampli- 
tude at the opposite side. With a 1500 points grid the nu- 
merical losses are only 5% for a wavelength of 1/30 of the 
computational box length , which correspond to 50 points 
per period. For larger wavelength the losses are smaller (2% 
at 1/20 and less than 1% at more than 1/10). As will be 
shown latter, spatial frequencies higher than 1/20 are not 
relevant to the problem considered here, thus the code is 
actually well suited for the present study. 

2.4. Nuclear Network and Detonations 

Since we are not concerned here with the details of nuclear 
burning, we adopt the classical a chain network, composed 
of 13 nuclei from ^^C to ^^Ni plus ^He. This time saving 
network is sufficient to describe the energetics of the C-l-0 
combustion (Timmcs 1999). For the reaction rates we use 
the NACRE compilation. The implementation of nuclear 
reactions in the hydro code HERACLES, uses a decoupled 
combustion step with multiple local nuclear time steps. To 
check this implementation we compare forced detonations 
with a high resolution adaptive mesh code, which can nicely 
resolve flames and detonations (ASTROLABE code, devel- 
oped by one of us: J. -P. Chieze). This is visible on Fig.(l), 
where the cells from ASTROLABE are plotted in blue dots 
along the Rayleigh line for mass and momentum conserva- 
tion, the Rankine-Hugoniot curve for energy conservation 
and the Crussard curve (the Crussard curve is the non- 
adiabatic Rankine-Hugoniot curve, modified by the energy 
released by the combustion). We can see, according to the 
ZND model, that the detonation results of a primary shock 
which heats up and compresses the plasma, pushing it on 
the Hugoniot. Then, combustion on scales much larger than 
the shock thickness, will expand the plasma and reach the 
intersection of the Crussard curve and the Rayleigh line. 
For the same set up, we plotted the simulation cells for the 
HERACLES code. Due to lower resolution, the initial state 
at V — 2 10~^ cm^ goes almost directly to the final 
state where the combustion is complete. But even without 
resolving the combustion zone, the final state corresponds 
nicely with the high resolution simulation . Then we com- 
pare several detonations from HERACLES and found deto- 
nation speeds agreeing within less than 1%. Finally we con- 
sider the effect of resolution on detonation initiation. The 
mechanism described later is self-consistent and rely only 
on compression and heating by shocks. There is no artifi- 
cial criteria or threshold for ignition. Usually only one cell 
runs away, when the density and temperature of the cell are 
sufficient so that heating due to combustion is faster than 
the cooling following the shock. In those conditions the cell 
reaches sufficiently high pressure to compress the neigh- 
bouring cells, propagating the detonation. This means that 
the ignition scale is not resolved. But since density and tem- 
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Fig. 1. Pressure (P)-specific volume iV) diagram of a 
forced detonation in a plasma of density po = 5.10^ gcm~^. 
The Hugoniot curve represents the energy conservation 
through a shock, the Rayleigh line is the line of mass 
and momentum conservation and the Crussard curve corre- 
sponds to the energy conservation after the nuclear burning 
is complete. Each blue point is the {P,V) value of a simula- 
tion ceU from the high resolution code, ASTROLABE. This 
simulation resolves the reaction zone following the leading 
shock, while the crosses representing the simulation cells 
of the code HERACLES, jump directly to the final state. 
This demonstrates that even without resolving the combus- 
tion zone, this code reproduces the correct energetics and 
dynamics. 



perature are averaged over the cell, the under-resolved igni- 
tion conditions are a lower estimate. Indeed, ignition occurs 
in the region of decompression following the shock, where 
temperature and density decrease almost linearly from the 
post shock state. Thus inside the detonating cell, the unre- 
solved region of maximum density and temperature would 
runaway faster than the whole cell, with averaged density 
and temperature. 



3. Shock formation and heating mechanisms 

In this section we analyse the propagation and the am- 
plification of sound waves in slab geometry. We consider 
a vertical column, with h measuring the height from the 
base {h = 0) to the top {h = L). The gravity points to- 
ward the base. By an analogy to spherical geometry, we 
speak there of inner boundary and upper density, while the 
top of the column would correspond to the outer boundary 
and lower densities. We adopt this simplifying assumption 
in order to derive a simple analytical description which is 
then compared to numerical simulations of hydrodynamics 
to validate them. 

3.1. Analytical description 

The Poynting flux of a wave of velocity amplitude uis F = 
^pu^Cg- This corresponds to the energy flux carried by the 
wave in the linear regime. This flux is conserved in the 
absence of dissipation. If the perturbation originates from 
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a region with density po and sound speed Cs,Oj its amplitude 
at any height h reads as: 



u{h) — uo< 



Pot's 



p{h)Cs{h) 



(2) 



In degenerate matter the sound speed depends mainly on 
the density. Thus a wave travelling outward will slow down, 
so that the perturbation amplitude will grow to conserve 
the flux. And this, until a shock forms and dissipates ki- 
netic energy. A good way of estimating the height of shock 
formation, Hgh, is to follow the hills of the wave until they 
catch up the valleys. The hills move with respect to the 
wave at -\-u{h) while the valleys at —u{h) and they will 
merge when they have travelled their initial separation of 
Ao/2. For a perturbation initiated at ft, = 0, this reads as: 



I poCsfi dh 



Ao 
4 ' 



(3) 



The height of shock formation is calculated for a baromet- 
ric structure with an upper density of po — 5.10^ g cm~^ 
decreasing to a density of 10"' g cm'^, over a length of 
L = 5260 km. The choice of the lower density is a limit 
where thermonuclear ignition would no longer be possible. 
The constant gravity acceleration is: g — GM©/ (1500A:to)^ 
and the sound speed at /i = is Cs^o — 10500 km s~^. 
Equation (3) for this gradient is plotted versus the fre- 
quency in Fig. 2. It provides a quite good estimate of the 
height at which a perturbation starts to dissipate and heat 
up the medium, due to the formation of shocks. There are 
two trends: the longer the wavelengths the farther they 
start dissipating and for a given wavelength Aq, stronger 
initial perturbation starts dissipating earlier. Hgh deter- 
mines if heating will occur in regions where densities are 
sufficient for ignition to occur, when Hgi^ < L. In that case 
the heating rate is increasing with the perturbation am- 
plitude u and the detonation initiation only necessitates a 
sufficient number of shocks to reach ignition temperature. 
On the other hand when Hgh > L, shocks form too late in 
regions of low densities, where detonation is not possible 
whatever the number of shocks. 

The amount of energy dissipated by a shock depends 
on its strength defined as 77 = (p* — p)/p, with p* and 
p the shocked and pre-shock densities respectively. To es- 
timate the shock strength we follow Ulmschncidcr (1967) 
and Ulmschncidcr (1970). Since we consider very weak per- 
turbations the resulting shocks are never strong {rj ~ 1), 
this implies they will be of the sawtooth type (Landau & 
Lifshitz (1959) §102 p385). For this shape the energy car- 
ried by the wave is: 

F=^ j\p - Po)udt = ^APAu, 

where AP and Au are the jumps in pressure and velocity 
across the shock front and T is the period. Assuming weak 
shocks, Au = rjCs and AP = pC^rj and using the adiabatic 
index 7 = 1 + P/E and = 7P/ p, this formula reduces 
to: 
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Fig. 2. Height of shock formation versus the initial wave- 
length of the perturbation normalised to the size of the box: 
Xo/L in the planar and spherical cases. 



On the other hand, shocks dissipate kinetic energy into in- 
ternal energy. The corresponding rate can be estimated fol- 
lowing Landau & Lifshitz (1959) for weak shocks, or for 
arbitrary strong shocks as in Brinklcy & Kirkwood (1947): 



1 dF 
Flh 



7 + 1. 



(5) 



where v — 1/T is the frequency. Then deriving eq.(4) and 
equating it to dissipation losses, we get the shock strength 
equation: 
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Osterbrock (1961) showed that for weak shock the con- 
served quantity is actually FC^. The inclusion of this re- 
fraction effect adds an additional factor ^ . Eq. (6) can 

then be integrated to get shock strengths 'q[h) given the 
position of formation and an initial shock strength 770. We 
note that the pressure term, ^g/C"^, and the energy dissipa- 
tion rate, (7 -I- l)rjv/Cs, dominate the RHS of equation (6), 
so that an asymptotic behaviour of the shock strength can 
be approximated by equating these two terms to obtain: 



19 



(7+ 1)^.7^ ■ 
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This relation is plotted as solid lines on Fig. 3 for different 
frequencies or equivalently initial wavelengths. It can be 
noticed that the asymptotic value depends only on the fre- 
quency of the perturbation. Accordingly perturbations of 
given wavelength with different amplitudes tend to a com- 
mon asymptotic behaviour: rjoc.uih). 

3.2. Numerical description and validation 

Here we use (6) to validate the numerical code used in the 
following section, which will be devoted to ignition. In our 
simulations, to create the initial sound waves we impose os- 
cillating density and momentum perturbations at the inner 
eulerian boundary, h = 0: 
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h (km) 

Fig. 3. Profiles of shock strengtfi. Dots represent tlie mea- 
surements from our numerical simulations: three different 
amplitudes (200, 100, 50 km s~^) in magenta, green and 
blue respectively for each wavelength families represented 
by their common asymptote eq.(7) for the three wave- 
lengths (Ao = 0.05L, O.IL, 0.2L in solid lines in red, ma- 
genta and black respectively). Also plotted arc the predic- 
tions from eq.(6) in black dashed lines. 
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Fig. 4. Propagation and amplification of a 100 km s^^ or 
M=0.01 perturbation in the planar density gradient. In blue 
the velocity from the numerical simulation, compared to the 
amplification of eq.(2) in the absence of dissipation. They 
separate when the shocks kinetic energy is dissipated. This 
is stressed by the temperature at three subsequent times 
plotted in black. 



( p{h = 0) = pq{1 - uosin{2Tri't)/Cs,Q), 

The velocity perturbation (uq) and the wavelength (Aq — 
Cs.o/v) are the free parameters of our study. For the outer 
boundary we constructed a transmitting boundary condi- 
tion, which allows matter carried by the shocks to leave 
the box. Fig. 3 summarizes our validation results and shows 
that the HERACLES code nicely reproduce the results of 
eq. (6) and (7) . These calculations have been performed 
for a set of 3 initial wavelength ranging from L/4 to L/20, 
where L is the length of the computational domain, each 
one with three different initial velocity amplitudes (200, 100 
and 50 km s~^). 

4. Ignition conditions 

The analytical results from the previous section can predict 
quite accurately the shock strength and the dissipation at 
any position in the degenerate carbon oxygen matter of a 
white dwarf. But, ignition by itself depends on the ther- 
monuclear properties of the plasma. Hereinafter we con- 
sider numerical simulation including the nuclear network 
described in section 2.4, both planar and spherical geome- 
try will be discuss. 

4.1. Slab geometry 

We showed that the perturbation wavelength plays a 
major role regarding the ignition process. Indeed, short 
wavelength perturbations generate a train of small am- 
plitude shocks whose cumulative dissipation brings the 
temperature up to ignition, at some point in the profile. 
Conversely a large wavelength perturbation can trigger ig- 
nition through dissipation across one single shock. To show 



this, we adopted an initial small Mach number perturba- 
tion, M=0.01, corresponding to a velocity amplitude of 
Mq = 100 km and we varied the wavelength by one 
order of magnitude. 

We consider first the short wavelength perturbation 
with Ao = 0.05-L = 263 fcm, corresponding to a frequency 
of 40 Hz. Transition to a shock occurs at the height of 
about 2500 km as predicted by eq.(2). This is the point 
where the velocity from the simulation, in blue on Fig. 4, 
diverges from the prediction of eq.(2) in red, since it enters 
the dissipation dominated regime. At that point the wave 
begins to heat up the plasma by dissipating kinetic energy. 
This is illustrated by the three black curves representing 
the temperature at three subsequent times. In that case, 
heating occurs through many successive weak shocks each 
one contributing to a small net temperature increase. This 
quasi-continuous dissipation process is governed by the evo- 
lution of the shock strength in eq.(5), nicely predicted by 
eq.(6) (see Fig. 3). For the a-chain thermonuclear network 
the runaway temperature is reached after 6.9 s (250 peri- 
ods). 

We turn now to the case of larger perturbation wave- 
length Aq — 0.5L, with the same initial amplitude, on Fig. 5. 
Detonation is triggered by dissipation of a single shock, 
before decompression occurs. The shock forms at about 
h — 4500 km, which again is close to the predictions of 
eq.(2). Due to a longer travel across the steep gradient, 
the perturbation is more amphfied, according to eq.(2), be- 
fore dissipation inhibits its growth. This leads to a much 
stronger shock, Au = 4000 kms~^, and thus a very strong 
dissipation. On Fig. 5, when the shock forms, the post shock 
temperature is greater than lO^K. But contrary to the 
short wavelength case it is immediately followed by drastic 
cooling in the decompression region. For ignition to occur, 
the post-shock induction time should be smaller than the 
cooling time behind the shock, which requires higher post 
shock temperatures. Since these single shocks are strong. 
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Fig. 5. Propagation and amplification of the same initial 
perturbation in the same gradient. Here the shock is formed 
latter and the numerical velocity (in blue) sticks to the dis- 
sipationless prediction (in red) much longer. This results in 
strong dissipation as soon as the shock forms and a deto- 
nation is triggered before decompression can occur. 
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Table 1. Parametric study of the parameter space \osc x 
with the conservative code HERACLES, for the case of a 
shallow gradient over 5000km with = GMQ/{1500km)^ 
and a steeper gradient with ten times the gravity. The ta- 
ble shows the number of shocks needed for ignition while 
dash symbols represent cases where ignition is impossible 
because shocks form at too low density. 
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Fig. 6. Ignition and propagation of a detonation in the 
large wavelength case (Aq = 0.5L and u = 100 km s~^). 
Temperature and pressure profiles in black and red respec- 
tively are plotted at four successive times. 



this regime is not predictable by eq.(6), but due to long 
wavelength it is more easily simulated. 

Whatever the mechanism to reach ignition temperature 
at the hot spot, the runaway creates a strong overpressure 
sufficient to propagate a double detonation both upward 
and downward. Fig. 6 illustrates such a successful ignition 
of a detonation. 

We analyse now the results of the parametric study. The 
number of shocks needed to ignite a detonation is presented 
in table 1, for 7 different frequencies and amplitudes varied 
from 200 kms^^ to 30 kms~^ , corresponding to Mach num- 
ber from 0.02 to 0.003. As discussed above increasing the 
initial amplitude, uq, has two opposite effects. It increases 
the acoustic energy flux, which depends on Mq, leading to 
stronger shocks but dissipation occurs earlier preventing 
further amplification. Regarding detonation ignition, the 



stronger the perturbation the faster the ignition. Also at 
some point, when decreasing the amplitude, H^h enters re- 
gion where the density is too low for ignition to be possible, 
this corresponds to the dash symbols in 1. 

The wavelengths of the perturbations has a drastic ef- 
fect as seen by comparison of the two cases described before. 
Although it has no effects on the energy flux {F=yu'^Cs 
), it has a strong effect on the height of shock formation: 
higher frequency or shorter wavelength perturbations will 
start dissipating earlier leading to less amplification and 
weaker shocks. The energy flux being independent of the 
frequency, the same energy will diluted and deposited over 
wider area. On the other hand it makes it possible for 
weaker perturbations to trigger a detonation given more 
time by reducing the shock formation height and bringing 
it back in denser region where ignition is possible. For ex- 
ample a 30 kms~^ initial perturbation, would not detonate 
with a wavelength of Aq = 0.5L, because the shocks would 
form too late, but a smaller wavelength makes ignition pos- 
sible after many weak shocks. 

Last, the density gradient has no effect in normalised 
units. The parameter determining the gradient is the accel- 
eration of gravity: g. If g is increased by a factor 10 (as is 
done in the lower part of table 1), the same density gradi- 
ent will shrink to 10 times smaller box. Taking accordingly 
smaller wavelength the shock strengths and height forma- 
tions are the same. The frequency is thus higher by a factor 
10, leading to a higher rates of dissipation. But in term of 
dissipation per shock, it is also the same. Thus ignition of 
a detonation needs almost the same number of shocks. The 
small difference is due to smaller region at maximum tem- 
perature just after the shock. It shortens the time where 
reactions can run away before decompression. 

4.2. Spherical Structure 

We apply this ignition mechanism to white dwarf struc- 
tures. In spherical geometry the flux is still conserved but 
since the surface of the spherical wave grows as the per- 
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Fig. 7. pCgT^ versus the initial perturbation position tq, 
the maximum of this curve is the position where a pertur- 
bation will be most amplified. 



turbation growth is weakened. The conservation equation 
now reads: 
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Fig. 8. Propagation of a 500 kms^^ perturbation of period 
36 ms initiated at a tq = 100 km. The period corresponds 
to a wavelength of Aq = 0.25Rwd at the optimal radius: 
ro = 350 km. Initially geometry dampens the wave but 
then the gradient steepens and it is amplified. We compare 
the theoretical formula (9) for a dissipationless case to the 
simulation. 



which gives for the velocity perturbation at any radii: 

y p[r)Cs(r) r 

Thus the amplification has to be stronger than the spherical 
damping rg/r. This introduces a new parameter: the posi- 
tion ro of the perturbation emitter. The dependence on tq 
of the amplification through eq.(9), is shown on Fig. 7. At 
small ro geometrical damping is strong and reduce the fi- 
nal amplification, whereas at large radius, the variations 
in density and sound speed are smaller, leading to weaker 
amplification. The maximum of p(ro)Cs(ro)r§ is the opti- 
mum position between these two opposite effects, and it lies 
around 350 km from the centre for the structure of the 1.43 
Mq cold C-f-O white dwarf used here. 

The first model of supernova progenitor chosen is a cold 
{T — 10^ K) self-gravitating sphere of carbon and oxygen 
{Xo = Xc = 0.5). Two different models are constructed 
with a total mass of 1.43 Mq, with the same fixed outer 
pressure {Pext — 1,7.10^^ ergcm~^) which correspond to a 
density of p '-^ 10^ g cm~^, but one with only carbon and 
oxygen and the other with a thin helium layer at the sur- 
face. The central density is pc = 7, 5.10^ gcm~^ and sound 
speed: Cg = 12000 km with a radius of Rwd = 1400 
km. We choose such a dense structure because it presents a 
very strong density gradient, best emphasizing the amplifi- 
cation mechanism. We will study more realistic structures 
of pre-expanded white dwarfs in next section. To initiate 
perturbation on this spherical model, the radius ro, of the 
oscillating pressure prescription, generating the initial per- 
turbation is chosen. This could correspond to the position 
of the turbulent flame which would create pressure pertur- 
bations. At the optimal radius of 350 km the density is 
2.10^5 cm"^ and the sound speed is Cg = 10000 km s~^. 

On Fig. 8, the result for an initial perturbation of 500 
kms~^ [M ~ 0.05) originating close to the centre (ro = 100 



km) is shown to emphasize the geometrical effect in the 
central region. Our simulation, in blue, is compared to the 
predictions of formula (9) in red. Initially where the curva- 
ture effect is strong, the damping is more important than 
the amplification due to the relatively fiat gradient in the 
central region, and the amplitude of the perturbation can- 
not grow. But then, the gradient steepens while the effect 
of growing wave surface becomes smaller and the waves are 
amplified. 

As in the planar case, we studied the possibility of det- 
onation ignition for initial perturbations of different wave- 
lengths and amplitudes. Here we consider the period of the 
waves, because the wavelength is varying with the posi- 
tion ro. The results are shown in table 2 and we observe 
the same trends. However, due to spherical damping the 
amplification is reduced, so that the formation of shocks 
strong enough to ignite a detonation necessitates larger 
amplitudes. In the end a quite large initial perturbation is 
needed: uq = 200 kms~^, corresponding to a Mach number 
of Mq = 0.02. We show in the following that if an helium 
layer is considered, the initial perturbation may be reduced 
to 30 km or Mq = 0.003. 

Helium Layer: The presence of an helium layer at the top of 
the C+0 core favours a detonation. This layer is small so as 
to keep the explosion undisturbed. According to Woosley 
& Kasen (2011), the presence of an helium shell above the 
C-l-0 core would lead to unobserved peculiar light curves 
and spectra if the shell is larger than ~ O.O5M0. Thus a 
very thin helium layer of Mhc = O.OIM0 is considered. 
Since Helium combustion is strongly exothermic this layer 
detonates with weaker shocks. Once ignited, when they 
reach the C+0 core these He-detonations are then strong 
enough to initiate C+O-detonations. The set-up here is dif- 
ferent from the double detonation scenario, where an accre- 
tion induced detonation is initiated at the base of the he- 
lium layer and where the transmission of the detonation to 



7 



by amplification of acoustic waves in type la supernovae 



C. Charignon and J.-P. Chieze: DDT 



■7/1-1 { IcTTI Q 

\ r\j 1 1 li J J 


1000 


500 


400 


300 


200 


100 


Tosc = 72 ms 


1 


2 


3 


4 






-'- osc — IILj 


2 


5 


7 


10 






T — 1 S -m s 

osc IlL^J 


4 


13 


17 


34 


67 




T — Q TO q 

osc <J 11 Ley 


7 


34 


78 


154 


438 




+Mhc = 0.01A'/q 




5u {km s^^) 


1000 


500 


250 


100 


50 


25 


Tosc = 72 ms 


1 


1 


1 


1 






Tosc = 36 ms 


2 


1 


1 


3 


5 


20 


Tosc = 18 ms 


4 


5 


6 


9 


19 


78 


Tosc = 10 ms 


7 


34 


35 


47 







Table 2. Parametric study of the parameter space Tosc x 
Uq with the conservative code HERACLES, for the case 
of a cold C+0 white dwarf (upper table) and 0.01 Mq 
helium layer (lower table). Perturbations are initiated at 
the optimal radius (rg = 350 km), corresponding to the 
maximum on Fig. 7. The period Tq^c = 72 ms corresponds 
to a wavelength of Aq = O.bRwD at the optimal radius of 
350 km. Numbers in the table correspond to the number of 
shocks needed for detonation ignition. 
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Fig. 9. Initial structures of the three pre-expanded white 
dwarfs. The solid lines represent the temperature profiles 
while the dashed lines are the density profiles of these white 
dwarfs. 



the core may need to wait for the He-detonation induced 
shocks to convergence inside the core and ignite the car- 
bon. Indeed helium either ignites above the core/envelope 
boundary, which leave time for the inward propagating det- 
onation to build up large enough pressure to directly trig- 
ger the C+O-detonation or ignition occurs at the transi- 
tion with sufficiently compressed and heated carbon so as 
to ignite an He-outward moving detonation and an inward 
C-f 0-detonation. And even if this fails under some condi- 
tions, the usually invoked shock convergence (Fink ct al. 
2010) could still trigger the C-|-0-detonation, with some 
additional delay. The results of the simulations with an he- 
lium layer are summarised in the lower part of table 2. Now, 
the ignition of detonation necessitates a smaller perturba- 
tion of 30 km s-^ or Mq = 0.003. 

Pre-expanded structures: However, observations of the SN 
la nucleosynthetic yield of intermediate mass elements re- 
quire that DDT should occur in a pre-expanded WD, typ- 
ically 1 to 1.5 seconds after the initial thermonuclear run- 
away. Ideally, the previous calculations should be applied to 
a dynamically expanding white dwarf, with the perturba- 
tion moving with the flame front. However, the HERACLES 
code is presently unable to follow dynamically the expand- 
ing structure. Thus we postpone calculations with dynamic 
expansion for further studies. Nevertheless, as a first ap- 
proach of the behaviour of the present mechanism in an 
expanding star, we consider three pre-expanded WD struc- 
tures with densities at the flame front of p// =9 10*, 3 10* 
and 10* g cm~^, respectively. They are obtained with the 
ASTROLABE code, using its moving mesh capabilities, to- 
gether with an ADR flame (Khokhlov 1995), to follow the 
whole star expansion, starting from ignition. These three 
structures used in the following are shown on Fig. 9, where 
the temperature and density profiles are displayed. Since 
the densities at the perturbations initiation positions are 
different in these three cases, we have to consider the Mach 
number of the perturbations. 

Due to expansion, density gradients steadily decrease, re- 
sulting in weaker and weaker amplification of perturba- 
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Table 3. Parametric study of the parameter space M x 
Tosc, for the three cases of expanded C-fO white dwarfs 
with an 0.01 Mq helium layer. Perturbations are initiated 
at the position of the flame. Tosc is the perturbation period 
in second, which varies from 0.2 to 0.02 s, while M is the 
Mach number, in the range 0.2 to 0.02. Numbers in the table 
give the time in second necessary for detonation ignition by 
the perturbation of period Tosc and Mach number M. When 
these numbers are followed by the letter c, it means ignition 
occurs in the carbon, otherwise ignition takes place in the 
helium layer. 



tions. Consequently, at later times, quite stronger perturba- 
tions are required to trigger sufficiently prompt DDT. Here 
we only consider detonation ignition occurring in less than 
1.5 seconds, thus excluding weak perturbations. Indeed, 
weak perturbations must be treated in a dynamical set up, 
allowing for repeated shock heating from the start. Table 3 
summarizes the results for the three pre-expanded white 
dwarfs. For those structures the sound velocities at the 
flame front are Cs = 8023, 6720 and 5880 km s~'^ respec- 
tively. To get an idea of the wavelengths of the perturba- 
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tions we consider, we give the period corresponding to a 
wavelength of half the star radius: T0.5 — 0.14, 0.21 and 
0.28 s respectively. 

We find again the same trends, longer period and 
smaller amplitude perturbations turn to shocks farther 
away, so that at some point shocks form in the outer en- 
velope, with too low density for ignition to be possible. 
On the opposite, shorter period and stronger perturbations 
start to dissipate closer to the flame, leading to ignition 
in the carbon-oxygen core rather than in the helium shell. 
These cases are marked with a letter c in the table 3. In 
this case, corresponding to high frequency perturbations, 
the energy is deposited mainly in the carbon-oxygen mix- 
ture and the presence of an helium layer is of no relevance. 
Only the strongest perturbations, reaching at least 10% of 
the pressure, can trigger a timely detonation. Moreover, 
due to the short distance between the flame, where pertur- 
bations are generated, and the locus of shock formation, 
where the energy is dissipated, the ignition doesn't really 
rely on amplification by the density gradient, since density 
does not significantly change over such a short distance. In 
fact, it relies mainly upon the natural evolution of any per- 
turbation to steepen into shocks, so that the ignition times 
are almost independent of the structure. However, since 
the energy is deposited close to the flame, the fact that the 
flame is actually moving, may prevent the temperature to 
reach ignition conditions. But in any case, the preheating 
of fuel due to shocks will have strong effects on the flame 
when it reaches this pre-heated region. Conversely long pe- 
riod perturbations with smaller amplitudes could trigger a 
detonation in the helium shell, due to a longer delay for 
shock formation. As the star expands, the minimum am- 
plitude to detonate the helium layer increases, due to less 
and less amplification through a flatter and flatter gradi- 
ent. From M — 0.003 in the cold and dense initial white 
dwarf, the required triggering pressure perturbation, bP jP ^ 
are respectively 0.02, 0.03 and finally 0.05. 

5. Discussion 

Finally we can estimate the energy fiux needed for a cold 
white dwarf with a thin helium layer to detonate through 
the amplification mechanism as: 

F„ = 0.5poC^,oU^ ^ 2.5 10^^ ergcmr"^ 

with pq = 2 10^ g cm"^ , Cs,o — 10000 km and a per- 
turbation amplitude of M=0.005 (u = 50 km s"^). The 
question is then to find how and where this flux could be 
produced. An estimate of the energy flux generated by the 
thermonuclear combustion during the initial deflagration 
phase at the density p^, is: 

Fnuc = Vfipenuc = 8.7 10^^ erg cm'"^ 

The laminar flame velocity used here is: Vfi — 75 kms^^ in 
agreement with Timmes & Woosley (1992) and the nuclear 
energy released is: e„„c = 5.8 10"'^'' erg g~^, corresponding 
to the combustion of the Xq = Xc = 0.5 mixture. The 
laminar deflagration speed is used as a lower limit since it 
should be greatly enhanced by turbulence. But even consid- 
ering this limit, only 0.3% of the generated energy injected 
into sound waves would be sufficient for a successful ig- 
nition in the cold initial white dwarf. In Table 4 the same 
comparison is made for the pre-expanded white dwarfs. For 
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9.3 10^' 


4.1 10^' 


2.5 10^' 


Fnuc {erg cm~'' s"^) 


1.6 10^^ 


2.2 10^^ 


2.3 10^' 



Table 4. Comparison of the acoustic flux needed for deto- 
nation ignition to the nuclear energy generation of the cor- 
responding laminar flame in the three pre-expanded white 
dwarfs studied previously. 



these structures the required energy represents a more and 
more substantial part of the available nuclear energy: 6%, 
19% and flnally almost all the energy generated by a lam- 
inar flame. Accordingly, only fast turbulent flames could 
generate, at late times, enough energy for this mechanism. 

As stressed in the previous section, a strong constraint 
on perturbations frequencies and amplitudes would be the 
time scale over which the transition to detonation has to oc- 
cur. In the context of the delayed detonation model of SNe 
la, the detonation has to be triggered after an initial phase 
of expansion driven by the energy released by the deflagra- 
tion. Thus, it should not be too fast at the early stage of 
flame propagation, so that an early detonation is avoided. 
But it should be fast enough, so that the detonation can 
be triggered before the expansion gets too large. In the 
classical local DDT, many studies agree to say the detona- 
tion has to occur when the flame reaches densities of about 
2 lO'' g cm~^ . In this paper the DDT mechanism is non 
locals with a shift in space and time between the deflagra- 
tion, producing the perturbations, and the position where 
they trigger a detonation. Thus the critical density, derived 
from usual DDT, is not directly relevant to this mecha- 
nism. Taking into account the delay due to wave propaga- 
tion and shock heating accumulation, we have considered 
white dwarfs where the flame is still at densities above the 
classical critical density, so that during the shock heating 
period, the star will further expand and detonation will 
be initiated when the flame encompasses regions with ap- 
proximately the critical density. Obviously a serious draw- 
back of this study is to consider a static flame position, 
in a structure devoid of previous shock formation. Indeed, 
shock heating should occur right from the beginning, with 
evolving characteristics as the flame travel across the white 
dwarf. The final heating will thus be a complex interplay be- 
tween the moving flame and the history of shocks formation 
and dissipation. Thus, the next step after this preliminary 
study, will be to take properly into account the expansion 
of the star, along with a moving perturbation generator. 

Finally, we found that there are two kinds of pertur- 
bations that could successfully trigger a non-local DDT. 
Either long period and small amplitude perturbations, cor- 
responding to large combustion scale, like Rayleigh- Taylor 
mushrooms. Those perturbations will heat up the outer he- 
lium layer, igniting an helium detonation, given the right 
shock heating history. Or, conversely, high frequency and 
large amplitude perturbations, corresponding to repeated 
strong pressure perturbations, which could be associated 
with small scale turbulent combustion. In that case it is 
unlikely that there will be enough time for shock heating 
to reach ignition temperature, since shocks will form not so 
far from the flame. But in case of detonation ignition fail- 
ure, the fuel would have been efficiently pre-heated, which 
would have a positive feedback on the flame. 
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In the classical DDT through the gradient mechanism 
of Khokhlov (1991b), pressure perturbation has to build up 
directly to detonation pressure, and this necessitate over 
10 km. In our case a failed build up to only a fraction 
of that pressure would be enough, if sufficiently repeated. 
Moreover, in an attempt to get direct transition to detona- 
tion, Poludnenko (private communication) obtained strong 
pressure pulses of about 30% of the unperturbed pressure, 
repeated every eddy turn-over time. This perturbation cor- 
responds to column 2 of table 3 and might detonate, given 
the right conditions. To conclude, we propose a promising 
possibility for the revival of failed local DDT, as non-local 
DDT induced by density gradient amplification. 
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